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BY 


E.  A.  SKELTON 


Summary 

V^~'~^lThe  cylindrical  pipe  Is  composed  of  a  number  of  concentric  layers  of 
elastic  solids  and  fluids.  It  is  excited  by  time-harmonic  mechanical 
point  forces  or  acoustic  sources.  The  exact  linear  equations  of  elasto- 
dynamics.  visco-dynamics  and  acoustics  are  solved  to  obtain,  via  the  finite 
element  method,  a  matrix  equation  which  connects  interface  displacements 
and  external  forces  in  the  spectral  domain.  Far-fieid  sound  radiation  and 
the  wavenumbers  of  free-waves  are  obtained  from  this  relation.  Numerical 
results  include  a  comparison  of  sound  radiation  levels  obtained  from  the 
exact  theory  and  a  shell  theory. 
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UST  OF  SYMBOLS 

Cartesian  coordinates 
cylindrical  coordinates 
spherical  coordinates 

radial,  circumferential  and  axial  displacements 

transforms  of  displacements 

axial  and  circumferential  harmonic  wavenumbers 

stresses  within  layer 

transforms  of  stresses 

pressure  in  fluid 

transform  of  pressure 

amplitude  of  point  source 

amplitude  of  point  force 

external  stress 

transform  of  external  stress 

location  of  point  source/force 

radian  frequency  of  vibration  (*2trf> 

density 

sound  velocity  of  acoustic  fluid 
Lame'  constants  of  elastic  layer 
dynamic  coefficient  of  viscosity  of  fluid 
acoustic  wavenumber  *  u/c 
longitudinal  wavenumber  *  u>(p/(  k+2p) ) 1/2 
torsional  wavenumber  *  v(p/p)l/a 
Y<ka-aa>  lm(y)  >  0 

V<kLa-aa)  Im(yL)  >  0 

V<*Ta-aa)  Imfry)  >  0 
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Inner  and  outer  radii  ot  layer 


INTRODUCTION 


The  mathematical  analysis  of  wave  propagation  in  and  sound  radiation 
from  a  single-layered  uniform  pipe  is  usually  based  on  a  shell  theory  to 
represent  the  pipe  vibrations  and  the  wave  equation  to  represent  the 
interior  and  exterior  acoustic  pressures.  Recent  work  using  shell  theory, 
of  relevance  to  the  problem  to  be  considered  here,  includes  that  of  Fuller 
A  Fahy  [1]  who  have  increased  understanding  of  the  physics  of  wave 
propagation  in  fluid-filled  pipes  through  their  discussions  of  wavenumber 
versus  frequency  plots  and  the  energy  distribution  of  free-waves  between 
fluid  and  shell,  and  that  of  James  (2.3)  who  has  shown  that  the  dispersion 
relation  helps  to  promote  a  qualitative  understanding  of  vibration  and  sound 
radiation.  Only  a  limited  amount  of  work  which  uses  the  exact  linear  theory 
of  elasticity  to  represent  the  vibrations  of  a  fluid-filled  pipe  has  been 
published:  in  particular.  Lester's  derivation  [4]  of  the  dispersion  relation 
of  thick  walled  pipes  is  central  to  the  approach  adopted  here. 

Whilst  there  is  a  considerable  amount  of  literature  available  on  wave 
propagation  in  and  sound  radiation  by  plane  layered  media,  the  literature  on 
layered  pipes  is  less  extensive.  Included  in  the  former  is  the  work  of 
Pestell  &  James  (51  and  Spicer  16]  who  use  the  exact  linear  equations  of 
elasto-dynamics.  visco-dynamics  and  acoustics.  Included  in  the  latter  is 
the  work  of  Yeh  A  Chen  17]  who  study  the  dynamics  of  two  coaxial  shells 
which  are  separated  by  a  viscous  fluid,  and  that  of  Avilova  &  Rybak  16] 
who  develop  the  dispersion  relation  for  a  thin  multi-layered  shell. 

In  this  Memorandum,  the  methods  of  Lester  14]  and  Pestell  A  James  [5] 
are  combined  to  give  a  matrix  relation  from  which  far-field  sound  radiation 
and  the  wavenumbers  of  free-waves  may  be  calculated  numerically  by  Fortran 
programs.  In  addition  to  being  a  valuable  tool  for  studying  free-wave 
propagation  in  and  sound  radiation  from  coaxially  layered  pipes,  the 
programs  are  useful  for  studies  to  determine  the  range  of  applicability 
of  shell  theories. 


2.  PROBLEM  FORMULATION 


(a)  General 


The  layered  system  shown  in  Figure  1  is  composed  of  an  arbitrary 
number  of  concentric  elements  each  of  which  satisfies  the  exact  linearised 
equations  of  acoustics,  elasto-dynamics  or  visco-dynamics.  It  Is  excited 
by  time-harmonic  point  sources  which  are  located  in  the  acoustic  fluids,  or  by 
time-harmonic  point  forces  applied  at  the  Interfaces  of  the  elastic  elements. 
The  time  factor  exp  (-tort)  is  omitted  from  all  equations. 


<b)  Fourier  Transforms 


It  Is  convenient  to  represent  the  variables  ss  Fourier  transforms. 
For  example,  the  displacement  components  are  represented  by 


ur(r.<p.z) 

m 

cos(mp) 

m 

ur(r.  n.  a) 

u^(r.^.z) 

*  ( 1/2ir)  X 

sin(mp) 

Up(r.  n.  a) 

uz(r.«p.z) 

n=0 

cos(n<p) 

uz(r.  n.  a) 

—  - 

exp(  iaz)  da 


(2.  1) 


where  the  spectral  displacement  components  are  given  by  the  inverse  transform 


ur(r.  n. a) 

2* 

cos  (ny>) 

u  r(r,<p.z) 

u^(r.  n.  a) 

=  (en/27r)  f 

sin(mp) 

u^(r.e.z) 

uz(r.  n.  a) 

0 

cos(n-p) 

uz(r.*.z) 

exp(-laz)dzd^  (2.2) 


The  Fourier  transforms  facilitate  the  solution  of  linear  differential 
equations  by  reducing  them  to  algebraic  equations. 


(c)  Construction  of  the  System  Matrix  Equation 


The  spectral  dynamic  stiffness  of  the  layers  and  the  interior  and 
exterior  media  are  matrices  which  relate  spectral  stresses  to  spectral 
displacements  at  the  Interfaces.  These  matrices 

(S(n.a)l.  lSe(n.a)]  and  [S*(n.a)] 

0*0  3*3  3*3 

are  obtained  for  the  acoustic,  viscous  and  elastic  media  in  Sections  3-7. 

The  assembly  of  the  element  matrices  to  form  the  system  dynamic 
stiffness  matrix  relation 


(SM(n.a)HU(n.a))  -  [EM(n.a)]  (2.3) 

M*M  M*1  M*1 

is  a  standard  finite  element  procedure  which  reflects  continuity  of 
displacement  and  equilibrium  of  stresses  at  the  interfaces.  Computer 
implementation  is  straightforward.  The  matrices  (U(n.a)]  and 
(EM(n.a))  are  column  vectors  of  interface  displacements  and  externally 
applied  interface  stresses,  respectively,  of  dimension  M*(3N+3)  for  a 
system  of  N  layers.  The  external  spectral  stresses  are  defined  to  be 
positive  when  acting  In  the  positive  direction  of  the  coordinate  axes. 

With  reference  to  Figures  2  and  3.  this  sign  convention  requires  changes 
of  sign  before  assembly  to  the  elements  of  (i)  the  last  three  rows  of  all 
matrices  (S(n.a)]  and  (il)  the  matrix  [S*(n.a)l.  An  example  of  an 
assembled  system  matrix  for  a  two-layered  system  Is  shown  in  Appendix  B. 


(d)  Free- Wave  Propagation 


In  the  absence  of  external  forces  or  sources,  the  system  of  homogeneous 
equations  (2.3)  has  a  non-trivlal  solution  if  and  only  if  det[SM(n.a)] 
vanishes.  For  given  values  of  n  and  of  w  there  will  be  real  and  complex 
values  of  a  at  which  the  determinant  vanishes.  Real  values  of  a  are  the 
wavenumbers  at  which  free-waves  propagate.  Complex  values  of  a  represent 
evanescent  waves,  whose  amplitudes  decrease  exponentially  with  distance. 
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Here.  the  real  branches  only  are  of  Interest,  and  plots  of  real  wavenumber 
a  against  frequency  for  selected  circumferential  harmonics  n  are  called 
wavenumber-frequency  or  dispersion  plots. 


(e)  Far-Fleld  Pressure 


In  the  special  case  of  the  exterior  medium  being  an  acoustic  fluid, 
the  far-field  pressure  in  spherical  polar  coordinates  is  found  from  the 
stationary  phase  approximation  to  be  [9] 

•  Ua(n.  k.  cos0)exp(-in7r/2)cos(n<p) 

p(R.  9.<p)  *  -iwpc  exp(ikR)/7rR  £  - 

n=o  sin©  Hn'(ka.$in6) 

where  Ut(n,  k.  cos0) .  the  radial  spectral  displacement  of  the  outermost 
interface  is  obtained  from  the  solution  of  the  matrix  equation  (2.3). 
evaluated  at  a*k.  cos0. 


3.  THE  ACOUSTIC  FLUID  LAYER 


Figure  2  shows  a  section  through  an  acoustic  fluid  layer  whose 
Inner  and  outer  boundaries.  r=a  and  r* b  respectively,  are  subject 
to  spectral  pressures  [p(n.a)]  which  produce  spectral  radial 
displacements  (ur(n.ar)]. 

The  pressure  in  the  acoustic  fluid  layer  satisfies  the  scalar 
Helmholtz  equation 


vap+kap  =  0 

whose  general  spectral  solution  Is 

p(r.n.a)  *  Ax(n.  a)  Jn(?r)  ♦  Aa(n,  o)  Yptyr) 
Evaluation  of  the  pressure  and  the  pressure-displacement  relation 

dp(r,n.a)/dr  =  put2  ur(r.n,a> 
at  the  boundaries  gives  the  matrix  relations 


p(b.n.a) 

Jn<rb) 

Yn(Yb) 

A1(n.a) 

p(a. n.  a) 

Jn(7a) 

YnOa) 

A2(n.a) 

ur(b.n.a) 

7 

m  ■■■— — 

Jn'<Tb) 

Yn'(yb) 

A1(n.a) 

ur(a.n.a) 

pu2 

Jn'(T*> 

Yn'(re> 

Aj in. a) 

(2.4) 


(3.  1) 


(3.2) 


(3.3) 


(3.4) 


(3.  5) 
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from  which  lA1(n.a).A2(n.dt)]T  may  be  eliminated.  This  gives,  after 
making  use  of  the  Wronskian  relation 


Jn(x)Yn'(x)  -  Jn'(x)YnCx)  =  2/ttx  . 


the  matrix  relation  connecting  spectral  pressures  and  displacements  at 
the  acoustic  fluid  layer  boundaries  as 


p(b.  n.  a) 

pwa 

p(a.  n.  a) 

yW 

Jn(7b)  Yn'fo'at-Yfjt'yb)  Jn'(7a>  -2/nvb 

2/irya  Yn(*ya)  Jn'tyb)  -Jn(ya)  Yn'(yb) 


ur(b.  n.  a) 
ur(a.  n.  a) 


where 


W  =  Jn'(«yb)Yn'(‘ya)  -  Jn'(-ya> Yn'(-yb> 

It  is  convenient  when  assembling  the  individual  layer  matrices  into  the 
system  matrix  to  write  equation  (3.6)  in  terms  of  the  fluid  stresses  as 

IS(n.  a)  ltu(n.  a)  1  =  (f(n.a)l 
6*6  6*1  6*1 


where  the  elements  of  [S(n.a)l  are  identically  zero  except  for 
Su  *  pua(Jn'(7a>Yn(yb>-Jn(rb)Yn'Cya>>/W 

S14  *  2pw2/iry  bW 

S4l  =  -2p<J2/7T7aW 

S44  *  pwa(Jn(7a)Yn'(7b)-Jn'(7b)Yn(7a))/W 


4.  THE  ELASTIC  LAYER 


Figure  3  shows  a  section  through  an  elastic  layer  whose  inner  and 
outer  boundaries.  r=a  and  r* b  respectively,  are  subject  to  spectral 
stresses  (r(n.a)l  which  produce  spectral  displacements  Cu(n.a)]. 

The  displacements  in  the  layer  satisfy  the  vector  equation  of 
motion  (10J 


(X+m)V(V.u)  +  MVau  *  paau/3la 
which  may  be  reduced  by  using  a  Lame'  substitution 

u  «  yf  +  vAjr,  ir  »  (O.o, Q)  ♦  va(o.o.h) 

to  the  three  scalar  Helmholtz  equations 

VaF+kLaF  ■  0.  VaQ+kjaQ  -  0.  VaH-»>kTaH  ■  0 


(3.  6) 


(3.  7) 


(4.  I) 


(4.  2) 


(4.  3) 
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Lester  [4]  has  shown  that  the  Fourier  transform  solutions  of 
equations  (4.  2-4.3)  together  with  the  transforms  of  the  stress-displacement 
relations 


^rr 

=  X(V.  u)+2p3ur/3r 

Tr<e 

s  fi(  3u^/3r-(  1/r)  u^>+(  1/r)  3ur/3v>) 

(4.4) 

rrz 

■  p(3ur/3z+3uz/3r) 

allow  the  construction  of  the  matrix  equation 

(S(n,  a) HiJ(n,  a)  1  =  £P(n.a)HR(n.  air^fiftn.a)]  =  fr(n.a)]  (4.5) 

6*6  6*1  6*6  8*6  6*1  6*1 

connecting  spectral  displacements  and  stresses  at  the  boundaries  of  the 
elastic  layer.  The  elements  of  the  matrices  £P(n.a)l  and  (R(n.a)l 
are  given  in  Appendix  A. 


5.  THE  VISCOUS  FLUID  LAYER 


Figure  3  Is  also  applicable  to  a  viscous  fluid  layer  whose  inner  and 
outer  boundaries.  r*a  and  r=b  respectively,  are  subject  to  spectral 
stresses  (r(n.at)]  which  produce  spectral  displacements  (u(n.a)]. 

The  linear  equations  of  motion  satisfied  by  the  particle  velocities 
and  the  fluid  pressure  within  the  viscous  layer  [111 

-vp  +  pvau  +  (p/3)v(v.u)  *  p3u/3t 

3p/3t  ♦  pV.  u  ■  0  (5.  1) 

dp/dt  *  ca3p/3t 
may  be  reduced  by  means  of  the  substitution 


U  *  VF  ♦  Va£.  jr  «  (0.0. Q)  +  Va(O.O.H)  (5.2) 

to  three  scalar  Helmholtz  equations,  together  with  an  expression  for  the 
fluid  pressure,  viz. 

VaF  ♦  lwa/(cJ-4IWM/3p)JF  -  0 

VaQ  +  ( ipu/fi)  O  »  0  (5.3) 

VaH  '  (ipw/p)!’  ■  0 

p  ■  *p/w)V.  u  (5.4) 
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The  stress-velocity  relations  in  the  viscous  fluid  layer  are 


rrr  =  -p-(2/t/3)V.u  +  2p8ur/0z 

rr<>  a  M(3u^/ar  -  (l/riu^  ♦(  1/r)  0ur/0<p)  (5.5) 

Trz  =  M(dur/dz  +  0uz/0r) 


from  which,  using  equation  (5.4).  p  may  be  eliminated  to  give  the 
stress-velocity  relations  for  the  viscous  fluid  layer  in  a  form  similar 
to  the  stress-displacement  relations  for  the  elastic  layer 


Trr 

Tr<p 

’’rz 

By  identifying 


=  (ic2p/<u-2M/3)  (V.  u)  +  2udur/dr 

ss  p(d u^/dr  -  (1/r)u^  +  (  1/r)  dur/8y>)  (5.6) 


=  fi(dur/dz  + 

0uz/3r) 

X  * 

(ic3p/(j  - 

2p/ 3) 

kT2  * 

ipw/p 

(5.  7) 

ii 

j e 

u2/(c1  - 

4ifa>A/3p) 

in  equations  (5.3)  and  (5.6)  it  is  evident  that  the  matrix  relation 
between  spectral  velocities  and  stresses  at  the  boundaries  of  the  viscous 
fluid  layer  is  of  the  same  form  as  that  between  stresses  and  displacements 
at  the  boundaries  of  the  elastic  layer.  It  is 


or 


[P(n.o)](R(n.  a)  l_1Iu(n.  a)J 
6*6  6*6  6*1 


-iw(P(n,a)HR(n.a)J-l(u(n,a)l 
6*6  6*6  6*1 


(r(n.a)l 

6*1 

[f(n.  o)l 
6*1 


(5.8) 


where,  subject  to  equations  (5.7).  the  elements  of  the  matrices  (P(n.a)l 
and  [R(n.a)l  are  given  in  Appendix  A. 


6.  THE  INTERIOR  MEDIUM 


(a)  Acoustic  Fluid 


Figure  2  shows  a  section  through  a  cylinder  of  acoustic  fluid  whose 
boundary.  r*b.  Is  subject  to  the  spectral  pressure  p(b.n.a)  which 
produces  the  spectral  radial  displacement  ur(b.n.o). 

The  pressure  In  the  fluid  satisfies  the  scalar  Helmholtz  equation 
(3.1).  whose  general  spectral  solution  equation  (3.2)  Is  subject  to  the 
condition  that  the  pressure  at  the  origin  remains  finite,  requiring  that 

A2(n.a)  ■  0. 
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Evaluation  ot  the  pressure  equation  (3.2)  and  pressure-displacement 
relation  equation  (3.3)  at  the  fluid  boundary  gives  the  relations 


p(b.n.a)  =  A3(n.  a)  Jn(-yb) 

ur(b.n,a)  =  (•y/pw2)  Ax(n.  a)  Jn'(vb) 

from  which  Aj^n.a)  may  be  eliminated  to  give  the  relation  connecting 
spectral  pressure  and  displacement  at  the  boundary 

p(b.n.a)  =  pu2(Jn(7b) /•yJn'('yb)]ur(b.  n.  a) 


(6.  1) 


(6.  2) 


It  may  be  convenient  when  assembling  this  element  into  the  system  matrix 
to  write  equation  (6.2)  in  terms  of  the  fluid  stresses  as 

(S*(n.  a)](u'(n.  a))  =  [f'(n.a)]  (6.3) 

3*3  3*1  3*1 

where  the  elements  of  the  matrix  (S'(n.a))  are  identically  zero  except  lor 

S*u  =  -p<J2Jn(7b) /yJn'(-yb) 


(b)  Elastic  Solid  or  Viscous  Muid 


Figure  3  shows  a  section  through  an  elastic  cylinder  whose  boundary. 
r=b.  is  subject  to  spectral  stress  [r*(n.a)]  which  produces  spectral 
displacement  ru'(n.a)l. 

The  displacement  in  the  elastic  cylinder  satisfies  the  vector  equation 
of  motion  (4.  1)  which  may  be  reduced,  by  means  ot  the  Lame  substitution 
equation  (4.2).  to  the  three  Helmholtz  equations  (4.3)  whose  spectral 
solutions,  after  imposing  the  requirement  that  the  displacement  at  the 
origin  remains  finite,  are 

F(r.n.a)  =  Ajfn.  o)  Jn(ri_r) 

Q(r.n.o)  =  A2(n,o)Jn(yTr)  (6.4) 

H(r,n.a)  =  A3  (n.  a)  Jn(y-rr) 

Evaluation  of  the  displacement  equation  (4.2)  and  stress-displacement 
relations  equations  (4.4)  at  the  cylinder  boundary  gives  the  matrix  equations 

Aj(n.a) 

lR*(n,a)J  Aj(n.a)  (6.5) 

3*3  Aj(n.a) 

A  x  ( n .  a) 

(r'(n, «)]  =  rP*(n, a)J  A2(n.at)  (6.6) 

3*1  3*3  A3 (n. a) 


Cu*(n,  a)) 
3*1 
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from  which  tA^n.a).  A2(n,a).  A3(n,at)]T  may  be  eliminated.  This  results 
in  the  matrix  equation  relating  spectral  stresses  and  displacements  at  the 
cylinder  boundary 

fSi(n.a)]lui(n.a)l  =  IP*(n.  a)  HR'(n,  a)  l-1[u*(n.  a) )  =  Ir'(n.a)]  (6.7) 

3*3  3*1  3*3  3*3  3*1  3*1 

where  the  elements  of  the  matrices  tPkn.a)]  and  lR*(n.  at)]  are  given 
in  Appendix  A. 

In  the  case  of  a  cylinder  of  viscous  fluid,  it  follows  from  Section  5 
that  the  matrix  equation  relating  spectral  stresses  and  displacements  at 
the  boundary  is 

-iw[P'(n.  at)  ](R'(n.  a)  l-1tu'(n.  a)]  =  [r*(n.a)l  (0.8) 

3*3  3*3  3*1  3*1 

where,  subject  to  equations  (5.7).  the  elements  of  the  matrices  (P'(n.a)l 
and  (R'(n.a)l  are  given  in  Appendix  A. 


7.  THE  EXTERIOR  MEDIUM 


(a)  Acoustic  Fluid 


Figure  2  shows  a  section  through  the  region  whose  boundary.  r=a.  is 
subject  to  the  spectral  pressure  p(a.n.a)  which  produces  the  spectral 
radial  displacement  ur(a.n.a). 

The  pressure  in  the  fluid  satisfies  the  scalar  Helmholtz  equation 
(3.1),  whose  general  spectral  solution  equation  (3.2)  is  subject  lo  the 
radiation  condition  that  at  large  values  of  r  the  pressure  consists  of 
outgoing  waves  only,  requiring  that 

Aa(n.  a)  =  iA^n.a) . 

Evaluation  of  the  pressure  equation  (3.2)  and  pressure  displacement  relation 
equation  (3.3)  at  the  fluid  boundary  gives  the  relations 


p(a.n.  a) 


A  jfn.  a)Hn(ya) 


ur(a.  n.  a) 


(y/pw2)  Ajtn.a)  Hn'(ya) 


(7.  1) 


from  which  A2(n.  at)  may  be  eliminated  to  give  the  relation  connecting 
spectral  pressure  and  displacement  at  the  boundary 


p(a.n.a)  =  pw2(Hn('ya> /yHn'(7a)  )ur(a,  n.  a) 


(7.  2) 


It  may  be  convenient  when  assembling  this  element  into  the  system  matrix 
to  write  equation  (7.2)  in  terms  of  the  fluid  stresses  as 
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(7.3) 


lS®(n,  a)Hu®(n,  a)]  =  lr®(n,a)l 
3*3  3*1  3*1 

where  the  elements  of  the  matrix  [S®(n,a)l  are  identically  zero  except 
for 


Sfx  =  -pw2Hn(ya) /yHn'tya) 


(7.  4) 


(b)  Elastic  Solid  or  Viscous  Fluid 


The  case  of  an  exterior  elastic  or  viscous  fluid  medium  is  not 
relevant  to  the  problem  of  far-field  radiated  sound,  but  it  is  included  for 
completeness.  Figure  3  shows  a  section  through  the  region  whose  boundary. 
r=a.  is  subject  to  spectral  stress  lf®(n.a)l  which  produces  spectral 
displacement  [u®(n.a)l. 

The  displacement  in  the  elastic  region  satisfies  the  vector  equation 
of  motion  (4.  1)  which  may  be  reduced,  by  means  of  the  Lame'  substitution 
equation  (4.2),  to  the  three  Helmholtz  equations  (4.3)  whose  spectral 
solutions  after  imposing  the  radiation  condition  of  outgoing  waves  only 
are 


F(r.n.ar)  ■  A2(n.  a)  Hn(?|/) 

G(r.n.a)  =  A2(n.a)Hn(yjr)  (7.5) 

H(r.n.a)  =  A3 (n,  a)  Hn(yrr> 

Evaluation  of  the  displacement  equation  (4.2)  and  stress-displacement 
relations  equations  (4.4)  at  the  boundary  results,  as  in  Section  6.  in 
the  matrix  equation  relating  spectral  stresses  and  displacements  at  the 
boundary 

CS®  ( n .  a)  J(u®(  n .  a)  1  =  [P9(n.  a)][Re(n.  a)  r1lu9(n.  a) )  =  {f®(n.a)I  (7.6) 

3*3  3*1  3*3  3*3  3*1  3*1 

where  the  elements  of  the  matrices  (P®(n,a)J  and  lR®(n.a)J  are  given  In 
Appendix  A. 

In  the  case  of  an  exterior  viscous  fluid  It  follows  from  Sections  5 
and  6  that  the  matrix  relation  between  spectral  stresses  and  displacements 
at  the  boundary  Is 

-lwlP®(n.o)HR®(n.a)]_1lu®(n.a)l  =  [T®(n,a)l  (7.7) 

3*3  3*3  3*1  3*1 

where,  subject  to  equations  (5.7).  the  elements  of  the  matrices  fP®(n.a)I 
and  [R®(n.a)l  are  given  In  Appendix  A. 


8.  POINT  FORCE  EXCITATION 


A  point  force  excitation,  of  magnitude  F0.  applied  in  either  the 
radial  or  axial  direction  at  cylindrical  coordinates  (r0.O.z0)  on  an 
elastic  layer  interface  may  be  represented  mathematically  as  the  external 
stress 


E(<p.z)  =  F00(z-z„)0(<>)/r0  (8.1) 

from  which  the  spectral  excitation  is  obtained  via  the  inverse  Fourier 
transform  equation  (2.2)  as 

E(n.a)  =  (enF0/27rr0)exp(-iaz0)  (8.2) 

The  addition  of  this  spectral  force  into  the  right-hand-side  of  the  system 
matrix  equation  (2.3)  is  a  straightforward  finite  element  procedure. 


9.  POINT  SOURCE  IN  ACOUSTIC  FLUID 


(a)  Green's  Function 


The  effect  of  a  point  source,  located  at  cylindrical  coordinates 
(r0.O.zQ)  in  the  acoustic  fluid,  is  best  developed  via  the  expansion  of 
its  free-space  Green's  function  (121 


p0exp(ikR0)/R0 


(ip0/2)^  encos(ny>)  J  Jn(yr0)  Hn(‘yr)exp(ia(z-z0)]da 
n=o  -• 


r>r0 
(9.  1) 


*  (ip0/2)£  encos(ny»  j  Jn(yr)  Hn(-yr0)exp[ia(z-z0)]dar 

n*o  “• 


r<r0 

(9.2) 


which  is  augmented  when  boundaries  are  present  by  a  scattering  term  which 
satisfies  the  homogeneous  Helmholtz  equation  (3.1).  Thus  the  total  pressure 
is 

p(r,«p.z)  *  p„exp(lkR0) /R0  +  ps(r.*,z)  (9.3) 

where  the  amplitude  of  the  scattered  wave  p8  is  determined  by  the  boundary 
conditions. 


(b)  The  Fluid  Layer 


The  spectral  form  of  equation  (9.3).  representing  a  source  in  an 
acoustic  fluid  layer.  Is  obtained  from  equations  (9.1),  ( 9.2 )  and  (3.2)  as 
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p(r.n.o)  =  (ip0/2)enJn(-yr0)Hn(yr)exp(-iatz0) 


♦  Ajtn.  a)  Jn(>r)  +  Aa(n,  a)  Yn(yr)  r*r0 

p(r.n.a)  =  (ip0/2)enJn(<yr)Hn(<yr0)exp(-iaz0) 

+  Aj(n.  a)  Jn(Yr)  +  A2(n.a)  Yn(-yr)  r<r0 

Evaluation  of  the  pressure  equation  (9.4)  or  0.5)  and  pressure- 
displacement  relation  equation  (3.3)  at  the  layer  boundaries  results  in  the 
matrix  equation  relating  spectral  displacements  and  stresses  at  the  layer 
boundaries 

IS(n.  a)][u(n.  a)l  =  I  Tin. a)}  +  tE(n.a)] 

6*6  6*1  6*1  6*1 

where  the  elements  of  (E(n.a)l  are  identically  zero  except  for 

Ej  *  2ip0enexp(-iaz0)CHn(7r0)Jn'(7a)-Jn(,yr0)Hn'(*ya))/W*Kb 

E4  *  2ipaenexp(-iaz0)(Hn(<yr0)Jn'(‘yb)-Jn(<yr0)Hn<Cyb))/W’ya 

Equations  (9.6)  and  (9./)  together  with  the  sign  convention  for  external 
stresses  show  that  the  point  source  in  the  acoustic  fluid  layer  has  an 
effect  on  the  right-hand-side  of  the  system  matrix  equation  equivalent 
to  external  radial  spectral  stresses  of  Ex  and  -E4  at  the  upper  and 
lower  boundaries  respectively. 


(c)  The  Interior  Fluid 


The  spectral  form  of  equation  (9.3).  representing  a  source  in  the 
interior  acoustic  fluid,  is  obtained  for  the  region  r>r„  from  equations 
(9. 1)  and  (6. 1)  as 

p(r.n.a)  =  ( ip0/2) enJn(yr„) Hn(7r) exp(-iazQ)  ♦  Aa(n, a) Jn(-yr) 

Evaluation  of  the  pressure  and  pressure-displacement  relation,  as  in 
Section  6.  results  in  the  matrix  equation  relating  spectral  displacements 
and  stresses  at  the  boundary  as 

(SUn.aOHu'fn.a)]  =  (rUn.a)]  ♦  lE'fn.a)) 

3*3  3*1  3*1  3*1 

where  the  elements  of  (E^n.a)]  are  identically  zero  except  for 

E[  *  2p0enexp(-iaz0)  Jn(7r0) /ybJn  (7b) 

Equations  (9.9)  and  (9.10)  show  that  the  point  source  in  the  Interior 
acoustic  fluid  has  an  effect  on  the  right-hand-side  of  the  ^system  matrix 
equation  equivalent  to  the  external  radial  spectral  stress  E\  at  the 
boundary. 


(9.4) 


(9.  5) 


(9.  6) 


(9.  7) 


(9.  8) 


(9.  9) 


(9. 10) 
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<d)  The  Exterior  Fluid 


The  spectral  form  of  equation  (9.3).  representing  a  source  in  the 
exterior  acoustic  fluid.  Is  obtained  for  the  region  r<r0  from  equations 
(9. 2)  and  (7.1)  as 

p(r.n.a)  =*  (ip0/2)enJn(<yr)Hn(‘yra)exp(-iaz„)  ♦  Ax(n,  a)  Hn(rr)  (9.11) 

Evaluation  of  the  pressure  and  pressure-displacement  relation,  as  in 
Section  7.  results  in  the  matrix  equation  relating  spectral  displacements 
and  stresses  at  the  boundary  as 

IS°(n,  a)  Hue(n.  a)  ]  *  [re(n.a)J  +  (Ee(n.ce)l  (9.12) 

3*3  3*1  3*1  3*1 

where  the  elements  of  IEe(n,a)]  are  identically  zero  except  for 

I®  =  -2p0enexp(-iaz0)Hn(>r0)/'yaHn'(>a)  (9.13) 

Equations  (9.  12)  and  (9.  13)  show  that  the  point  source  in  the  exterior 
acoustic  fluid  has  an  effect  on  the  right-hand-side  of  the  system  matrix 
equation  equivalent  to  external  spectral  radial  stress  of  -E^  at  the 
boundary. 


10.  NUMERICAL  RESULTS 


(a)  General 


Fortran  programs  have  been  written  to  calculate  and  plot  the  axial 
wavenumbers  of  free-vibration  and  the  far-fleld  sound  level  due  to  point 
source  or  point  force  excitation.  The  programs  are  written  in  double 
precision  complex  arithmetic  which  Is  simulated  [6]  by  the  use  of  double  precision 
variables  of  leading  dimension  2.  They  were  run  on  a  PDP-I1/34A  computer 
for  which  the  arithmetic  word  length  of  32  bits  used  for  single  precision 
complex  arithmetic  Is  inadequate  in  this  type  of  problem. 

The  material  and  geometric  constants  In  SI  units  which  are  used  in  the 
calculations,  involving  a  single  pipe  and  two  coaxial  pipes  separated  by 
water,  are  as  follows: 


Steel  inner  Pipe: 

X=10. 44  x 
a*0. 2096 

10l° 

b»0. 

p-7. 

2350 

56  x 

10l° 

p-7700.0 

Steel  Outer  Pipe: 

X*10. 44  x 
a*0. 2700 

10l° 

b»0. 

p«7. 

3000 

56  x 

10l° 

p-7700. 0 

Water: 

p-1000.0 

C-1500.0 

Air: 

p*l . 21 

c*343.  0 

Damping  In  the  pipe  walls  Is  included  by  setting  X  and  u  to  the  complex 
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values  XO-ii?*)  and  fid- ii)^) .  respectively,  where  the  numerical  values 
of  the  hysteretlc  loss  factors.  Vx  and  are  chosen  as  0.02. 


(b)  Wavenumber  Plots 


Figure  4  shows  the  real  branches  of  the  wavenumber  versus  frequency 
plots,  for  the  first  three  circumferential  harmonics,  of  the  inner  pipe 
containing  water  and  surrounded  by  a  vacuum.  The  plots  are  identical  to 
those  obtained  by  Lester  (41  who  has  also  superimposed  plots  obtained  from 
a  shell  theory.  The  physical  interpretation  of  the  plots  is  given 
elsewhere  [11. 

Figure  5  shows  the  plots  for  the  composite  pipe  consisting  of  the  inner 
and  outer  pipes  separated  by  a  3.  50cm  layer  of  water;  again  the  interior 
is  water  and  the  exterior  a  vacuum.  The  plots  contain  many  more  branches 
than  the  previous  plots  because  of  the  partial  Independence  of  waves  in  the 
solids  and  fluids. 


(c)  Sound  Radiation  Spectra 


Figures  6-11  show  the  variation  of  far-field  'airborne'  sound  pressures 
<dB  ref.  1  micropascal  at  1m)  with  frequency.  The  plots  (all  at  <p=0a) 
consist  of  straight  line  joins  which  connect  levels  computed  at  15Hz  spacing; 
hence,  the  sharp  peaks  are  unlikely  to  have  been  resolved  fully.  The 
excitation  is  either  a  IN  rms  radial  force  applied  to  the  inner  pipe  at 
the  cylindrical  coordinates  (0.2096.0.0).  or  It  is  an  interior  point  source, 
located  at  (0.1482.0.0),  whose  free-field  pressure  Is  120dB. 

(d)  Radiation  from  Inner  Pipe 


Figure  6a  shows  the  airborne  sound  radiation  from  the  inner  pipe  at 
6=90 •  due  to  point  force  excitation.  Figures  7a.  8a  and  9a  show  the 
sound  radiation,  at  8=90*.  80*  and  70*.  respectively,  due  to  point-source 
excitation.  The  corresponding  (b)  Figures  contain  the  sound  radiation 
calculated  from  a  shell  theory.  There  is  good  agreement  between  the  (a) 
and  (b)  spectra  up  to  the  ‘ring  frequency'  of  3.6kHz.  the  differences 
thereafter  being  mainly  attributable  to  the  relatively  Increasing  density 
of  the  peaks  In  shell  theory. 

James  (2)  has  demonstrated,  using  shell  theory,  that  the  frequencies 
at  which  'peaks'  appear  in  the  sound  spectra  depend  on  the  location  (9) 
of  the  observation  point.  These  frequencies  are  the  frequencies  at  which 
the  lines  a*k.cos(0)  cross  the  wavenumber  branches,  where  k  is  the 
wavenumber  in  air.  The  same  Is  true  when  exact  linear  theory  Is  used. 

In  fact,  the  difference  between  the  (a)  and  (b)  plots  Is  reflected  in  the 
differences  between  the  dispersion  plots  of  the  shell  and  exact  theories. 
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(«)  Radiation  from  Composite  Pipe 


Figure  10  shows  the  airborne  sound  radiated  at  0=90*  for  the  cases 
of  (a)  point  force  excitation:  and  (b)  point  source  excitation.  Figure  11 
shows  the  sound  radiation  due  to  point-source  excitation  calculated  at 
(a)  0=80*  and  (b)  0=70*. 

Again,  the  frequency  at  which  'peaks'  appear  in  the  spectra  may  be 
found  from  the  dispersion  plots  as  the  frequencies  at  which  the  lines 
a=k.cos<0)  cross  the  wavenumber  branches.  The  dispersion  relation  is 
equally  useful  in  interpreting  the  spectra  from  layered  pipes. 


E  A  Skelton  (SO) 
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APPENDIX  A 


The  Elastic  Layer  Matrices  tP(n.a>l  and  [R(n.a)) 


11 

Pn- 


IS 
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.  5. 


xi 


1  s 

R11 

R,3 

R». 

R_. 


l  f 


*.  5. 


6*6 

6*6 

2*rlJZ(rLb)-x*lJn(r Lb) 

Pn- 

2^A(V»>-K*£V1'Lb> 

( 2#n/b2 )  [ y (‘yTto)-JnC>T^>)  J 

pi." 

(2<m/b2)  [rTt>^CyTb)-YnO-,IJb)  ] 

-Ziainr'j-Cy^) 

P  - 

It 

-2iaio-2Y*(rTb) 

( 2nn/b  * ) [ Jn ( rLb ) -i^bj^ ( y^b ) ] 

P„‘ 

( 2n<i/b2 )  1 *n  ( >Lb )  -rLbY^  ( >Lb )  ] 

Oi/b2)[2-yTbJ|;(-yTb) 

Pi«- 

(M/b2)[2-yTbY|^CyTb) 

+  (-»^b*-2n,)Jn(‘yTb)} 

•*-(>^>2-2n2)YnCyTb)] 

( 2  i«<*n/b 2  w ) -Jn  ( yTb)  ] 

pa.’ 

( 2  ia/in/b  2 )  [brTY^  ( -y^ )  - *n  ( yjo  )J 

2i«MVAC’Lb) 

P„- 

2taa>Ly;(yLb) 

(  ian/n/b)  Jn(‘y<rb) 

P*4- 

(ian<i/b)Yn(-yTb) 

M-rT(2a2-k^)Ji;cyTb) 

P3.“ 

ai'T(2a2-)^,)Yi;(yTb) 

and  6  are  obtained  by  setting 

b=a 

in  rows  1.  2.  and  3. 

VA(>Lb) 

R13- 

VA(V» 

(n/b)Jn('yTb) 

Rl«- 

(n/b)Yn(rTb) 

“iaVA(V,) 

Rl.- 

-laVA(V,) 

-(n/b)Jn(7Lb) 

R»" 

-(n/b)Yn(TLb) 

“VA<V» 

R*4- 

“WV** 

( ian/b)Jn(7Tb) 

R»- 

( ian/b)Yn(yTb) 

iaJn(1,Lb) 

R)S- 

iaYn(rLb) 

0 

R,«" 

0 

R3.- 

and  6  are  obtained  by  setting 

b=a 

in  rows  1.  2.  and  3 
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The  Elastic  Interior  Matrices  IP^n.a)]  and  (R^n.a)] 

3*3  3*3 


"  2^Lj;(,'Lb)-XkLJn(l'Lb> 

-  (2^/b2)[7TbJ;i(7Tb)-Jn(7Tb)] 

-  -2iaMT^(7Tb) 

-  (2nM/b2)[JnCyLb)-7LbJ^(yLb)] 

-  (M/b2)[27TbJ^(7Tb)+(7jb2-2na)Jn(7Tb)] 

-  (2ia/m/b2)[b7TJ^(7Tb)-Jn(7Tb)] 

-  2i«M7LJi;(7Lb) 

-  (ia/in/b)Jn(7Tb) 

-  M>T(2a2-l4)J^(7irb) 


*LJA<*Lb> 

(n/bJJ^.jb) 

-ion,TJA(7Tb) 

-(n/b)Jn(vLb) 

"rTJn(rTb) 

(ian/b)Jn(7Tb) 

iaJn<TLb) 


(a2-k^)Jn(7Tb) 


♦ 


I 

I 


The  Elastic  Exterior  Matrices  lPe(n.a)J  and  (Re(n.a>) 
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APPENDIX  B 
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